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ABSTRACT 

Momentum feedback schemes, such as the one currently proposed 
for the Apollo Telescope Mount cluster, provide a relatively simple 
and effective means of obtaining near-linear control from highly non- 
linear control moment gyro (CM3) momentum storage systems. However, 
these schemes are subject to adverse orientations of the CMG momentum 
or spin vectors for which spacecraft control capability is lost, 
although the momentum storage capacity of the CMG system has not yet 
been fully used or saturated. The adverse CMG orientations occur when 
all CMG momenta and the desired command momentum are collinear and can 
be classed as saturation, anti saturation, parallelism, and 
antiparallelism. 

Since the adverse CM3 orientations generally occur only for 
rare combinations of the command momentum and the CMG orientations, a 
simple correction technique for these conditions appears desirable in 
order to minimize the onboard control computer requirements. A study 
was made, leading to development and analysis of a technique allowing 
the CMG momentum vectors to avoid the adverse orientations and thus 
to provide optimal control. The technique was formulated in terms of 
a feedback scheme using a check factor K formed as the product of 
the dot products of the three CMG momentum unit vectors. These are 
already available, as direction cosines, in the momentum feedback loop 
of the control circuit. The parameter K approaches unity for 
adverse orientations and is used to initiate a single correction 
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maneuver to rotate the CMS momenta away from a pending adverse 
orientation and to reduce K to an acceptable value. A digital 
computer program for AIM stations was modified to incorporate this 
correction maneuver and perfonoance was evaluated in simulated cases. 
Antiparallelism is the most typical adverse orientation expected and 
was used as the primary example for the evaluation of the feedback 
optimization technique. The effectiveness of the technique for other 
adverse conditions was also examined. 

The K-factor correction technique described in this paper 
permits avoidance of adverse CMG orientations up to saturation and 
produces control for exit from the adverse orientations in the event 
that they are arbitrarily acquired. 
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CHAPTER I 


INTRODUCTION 


The subject area is optimal feedback control for an advanced 
control system concept for long-term spacecraft. The control con- 
figuration is a momentum storage system using control moment gyros 
(CMG’s). A full-scale CMC unit is shown in Figure 1. A CM3 consists 
of a flywheel spinning at a constant speed and mounted on double 
gimbals schematically shown in Figure 1(a). 


H - Vector spin axis 


Torque r motor- 


inner gimbal axis 



Torquer motor 


Outer gimbal axis 


Flywheel 


Figure l(a).- CMC schematic. 

These gimbals are alined with a spacecraft axis in their 
reference position. The outer gimbal cannot move into gimbal lock 
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and rotates ±l8o° from a reference angle to take advantage of the 
available momentum storage capacity. The inner gimbal can incur 

grnbal lock at ±90° and is limited to an angle smaller than this 
value by means of mechanical stops. 

Control torques are introduced by the precession of the constant 
speed flywheel which is mounted on the gimbals. For example, if a 
torque is applied to the outer gimbal of the CMJ, an equal and opposite 

torque will act on the spacecraft, and the flywheel rotor will precess 
about the inner gimbal. 

The CM! has been shown to be an effective means of momentum 
storage for the stabilization of long-term manned spacecraft 
(Ref. 1). The first manned spacecraft to use CkG's for attitude 
control is the Apollo Telescope Mount (ATM) cluster configuration 
comprised of a Saturn S-IVB Workshop and a multiple docking adapter 
to which is attached the Apollo Command Service Module (CSM), and the 
Lunar Module - ATM rack which houses a complement of solar pointing 
experiments (Fig. 2). The control system utilizes three CM3's mounted 
orthogonally on the exterior of the ATM rack (Fig. 2(a)). Since 
maneuvers can best be carried out with near-spherical momentum 
envelopes, all three gyros have been selected to have the same 
angular momentum of H = 2000 ft-lb-sec. The requirement for the 
CMJ control system is to point the ATM cluster to at least 0.1° about 

all axes while a vernier platform points the experiment package to 
2*5 seconds of arc. 
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Figure 2.- Apollo telescope mount cluster configurati 
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The CMG system inherently exhibits nonlinearities in its 
response and requires angular momentum feedback to establish satis- 
factory linear control (Appendix A) . With the three distinct CMG angu- 
lar momentum vectors to oppose the commanded momentum, the momentum 
feedback scheme suffers from variable gain as a function of the relative 
orientation of the CMG momenta and the commanded momentum and is also 
subject to adverse CMG system orientations which lead to loss of con- 
trol and large spacecraft errors. Because of these f’MG system control 
characteristics, techniques are needed for minimizing system gain 
variation and avoiding adverse CMG orientations. A primary objective 
is minimization of inner gimbal angles to avoid the selected 80 ° limit 
and to reduce gimbal power requirements. 

One technique defined as the isogonal correction, employs a 
distribution law and rotation law for the individual angular momentum 
vectors (Ref. 2). The desirable distribution is for the individual 
angular momentum vectors to have equal components along the total con- 
trol momentum, and thus maintain the angular momentum vectors in an 
umbrella-like configuration which effectively prohibits the attainment 
of an adverse orientation. The rotation law is added to minimize the 
inner gimbal angles and avoid the inner gimbal limits as much as 
possible. The mechanization of this type of law requires significant 
additional onboard computer capacity to carryout the additional calcu- 
lations. The use of this complex procedure and added control computer 
capability appears undesirable for a flight mission since the adverse CMG 
orientations generally occur only' for rare combinations of the command 
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momentum and the CMC initial orientations. It was therefor* decided 

that a study of a simpler alternate method of first detecting and then 

correcting for the approaching adverse orientations would be extremely 
useful. 

The problem areas in the developnent of this simplified technique 
included: 

1* ^ determination of a check factor "K" as an indicator of 

momentum orientation. 

2. A deadband technique to keep the check factor within a 
given range to avoid the adverse orientations in order to keep the 
spacecraft errors within minimum acceptable levels. 

3. A method for varying K to achieve the approximate isogonal 
correction. 

These problem areas constitute the purpose of this research 
and are examined herein. A performance evaluation with digital 

computer simulations has been made with the best selected feedback 
technique to both avoid and exit from the adverse CM! orientations. 


CHAPTER II 


SPACECRAFT AND CONTROL SYSTEM SIMULATINN 

For ATM mission simulations a digital computer solves the 
spacecraft rigid tody and flexible equations of motion. The equations 
comprising this computer program are from Reference J and are included 
in the Appendices. The contribution of the flexibility to the space- 
craft rate and attitude errors is determined from dynamic and quasi- 
static modal representations of the AIM cluster included in Appendix B. 
For the present study only biased disturbance torques were applied to 
the spacecraft. From these disturbances and the resulting spacecraft 
rate and attitude errors, the computer then solves the control law to 
command the appropriate gimbal rates to the CMG torque motors which 
precess the spinning rotors and apply restoring control torques to 
close the simulation loop. Reference axes for the ATM cluster config- 
uration are defined in Figure 5 where the spacecraft 2 -axis is taken 

to be the solar pointing axis, and the x-axis to have the minimum 
moment of inertia. 

The control system for the spacecraft uses a gimbal rate 
steering matrix and CMG system momentum feedback as presented in 

Appendix A. The computer program constants used for the ATM simulations 
are included in Appendix C. 


POINTING AXIS 



ROLL 



Figure 3»- ATM axes system. 



CHAPTER III 


ADVERSE CMS ORIENTATIONS 

Digital computer simulations with CMS control systems have 
indicated serious control problem areas associated with particular 
orientations of the CMS angular momentum or spin vectors that could 
occur during normal control operations. 

The CMS axis placement for the ATM orientation and positive 

gimbal angle notation for the CMS angular momentum vectors is defined 
in Figure 4. 


z A 



X 


Figure 4. - CMS axes system. 
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variables 3j_ represent the CMi girnbai angles with the 

subscript identifying the CM3. The outer global angles are designated 
by the symbol a and the inner gimbal angles by 3 . The positive 

outer gimbal rate is about the -z spacecraft axis, about the 

-x spacecraft axis, and about the -y axis. 

Each of the three CMG units has an angular momentum of 


H = 2000 ft-lb-sec 

and limits are imposed on the gimbal rates and angles for the ATI.' 
system. The CMG inner gimbal angles aie limited to ±80° as previously 
mentioned, and the CM3 outer gimbal angles are limited to ± 175 ° with 
respect to this zero momentum orientation for which 


“l = ^5° 3 X = 0° 

= ^ 5 ° 3 2 = °° 

a 3 = ^5° 3 5 = o° 


The CM3 gimbal rate magnitudes 0 ± are restricted to be equai to 

or less than 3-5 deg/sec by gimbal power limitations. 

For the three-gyro system, there are four possible CMG spin- 
vector orientations which should be avoided durir* normal operations. 
Figure 5 indicates the relative placement of the disturbance moment, 

* and the CM:; an £ ular momentum vectors for the undesirable 
orientations. For these relative placements of the disturbance and 
control vectors, the CMG control system cannot produce a control torque 


* 
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(a) Antiparallelism 


(b) Saturation 



(c) Parallelism 


(d) Antisaturation 


Figure 5 .- Adverse CM3 spin vector - cotnnand moment orientations. 
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to oppose the disturbance torque and excessive spacecraft errors will 
occur. 

For the antiparallel condition bias disturbance torques have 
driven the CM3 spin vectors into a configuration with two H-vectors 
opposing the third and the disturbance moment, {%}, collinear, but 
opposed to the third CM3 spin vector. For the saturation condition a 
bias disturbance torque has driven all three CM3 H-vectors parallel to 
each other and alined with and having the same sense as the disturbance 
moment . The parallel condition is identical with the anti- 

parallel condition in that two H-vectors oppose a third, but the 
disturbance moment is now alined with the third CMC spin vector. 

For the antisaturation condition all three CM3 H-vectors are parallel 
to each other and oppose the disturbance moment. With CM3 control 
lost in any of these adverse orientations, the spacecraft mission is 
seriously handicapped, and a means of correction for the adverse 
conditions is needed. 

The orientations that predominate during normal operations are 
the saturation and antiparallel conditions. For the saturation con- 
dition all of the control system momentum has been expended and 
external moments frcrn either reaction jets or spacecraft reorientations, 
prior to total saturation, are required to change the CM3 spin-vector 
alinements with the bias command moment. In normal earth-orbit opera- 
tions, the gravity gradient bias moments will eventually saturate the 
control system, and planned desaturation techniques will be applied 
during the dark side or night portion of the 


A 


1 




orbit so as not to interfere 


with mission requirements. The antiparallel condition, on the other 
hand, has used only one-third of the control system momentum, even 
though effective control is no longer possible. 

For these pending adverse orientations, the effective gain of 
the control system is reduced since the control torque decreases as 
the adverse orientation is approached. Optimal feedback control can 
only be obtained by forcing the CMS vectors to maintain favorable 
orientations relative to the commanded momentum as long as possible. 
The associated fundamental momentum relations for the CM3 system are 
developed in the following sections. 


CM} Momenta 


Unit vectors for the individual CM3 angular momentum vectors 




-< h 3> 


where the vectors fa), fa), (hj) ccmprise the direction 
cosines for the CM} momenta and became 
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-cos P 2 sin ag 
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* 25 J 



-cos Pj sin a 
( -sin 
cos 3 ^ cos 
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( 2 ) 


The components for the h-y are available directly from the CM 3 
gimbal resolvers. 


Normalized Momentum 

The total control system angular momentum is 

( H ) = ( H l} + (Hg) + {** 3 } (3) 

or 

{ H )- H [{ h i) + ( h a> + (»3>] 00 

and a normalized total momentum magnitude is defined as 


HT = 


[(h) • (H)T/ g 

JH 

or 

m - i [( 0 *> ♦ to - to ) ■ (to ♦ to ♦ („,))] 1/2 

with the resulting bounds 


0 < HT < 1 


For the antiparallel and parallel orientations 


HT = 


1 

3 


and for the saturation and 


an oi saturation orientations 


HT = 1 


Antiparallel Orientation 

consider a pending antiparallel orientation examine the CM3 
spin-vector alinements sketched below. 


{%} 
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Two of the CMS spin vectors, here chosen to be (h 2 ) and (Hj) 
approach each other, while the third, ( k ± ) is opposed to the commanded 
momentum (h c ) . The antiparallel condition is reached when 

( H l ) • fe ) = -« 2 

{%} • (H 3 > = -H 2 ( 8) 

(h 2 ) • (h 5 > = H 2 

or in terms of the unit vectors when 

( h l) ' ' 'I 

(hj) • (hj) = -i (9) 

• (h 5 ) = 1 

Having developed an expression for the normalized momentum, HT, the 
next step is the determination of a check factor, K, as a nomadized 
indicator of momentum orientation and relate both normalized quantities 
to the adverse orientations. 



CHAPTER IV 


DETECTION OF ADVERSE ORIENTATIONS 

To detect an approaching adverse orientation the dot products of 
the direction cosine vectors are computed as 

k 12 = (hi) ' (h 2 ) = h u h 21 + h 12 h 22 + h 13 h 23 

K 1J ( h l) (hj) = hlihji + hi 2 h 52 + h-L^hjj ( 10 ) 

k 23 = ( h s) * (h 5 ) = h 21 h 51 + h 22 h 52 + h 23 h 53 

and a variable check factor K is defined as 

K = K 12 K 13 K 23 (11) 

From the reasoning leading to equation (ll), it is apparent 
that an adverse condition is pending if 


K ->l 


( 12 ) 


The normalized moment™ ffr and the check factor K 

the specific adverse orientations end hence establish 
control regions. 


together define 
undesirable 




* Dezween 




momentum HT for the four adverse condi 


K and the normalized 
is summarized in Table I. 
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TABLE I. - CHECK FACTOR - NORMALIZED MOMENTUM RELATION ^ 

FOR ADVERSE CMC ORIENTATIONS 


Adverse condition 

Check factor, K 

Normalized 
momentum, HT 

Antiparallelism 

1 

1/3 

Parallelism 

1 

1/3 

Saturation 

1 

l 

Antisaturation 

1 

l 


Linear Analysis of Impending Adverse Orientations 
With a constant disturbance torque bias (characteristic of 
external torques in an orbit environment) applied to the spacecraft 
the CMC control system mist eventually reach one of the collinear 
orientations. The objective of the check factor K reduction is to 
keep the control system away from pending antiparallelism and 
parallelism at HT = l/j, for the constant disturbance torque, and 
allow the system orientation to proceed to saturation at HT = 1 
*° T ’ e P lanned desaturation techniques will be applied. In this way, 
the full capability of the th. - louble-gimbaled CMS system may be 
realized prior to any adverse control effects on the spacecraft 
response. 

It was previously indicated that K = 1 for all adverse CM3 
orientations. The determination of how the total system momentum 
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varies with the position of the CMG spin vectors, especially for a 
pending adverse condition, is the necessary first step of analysis 

since the check factor K has been chosen as an indicator of momentum 
orientation. 

In correcting for a pending antiparallel orientation where the 
check factor K approaches 1.0, a deadband technique is analyzed and 
developed to reduce X. This requires the definition of a check factor 
limit, K^r^, such tnat when 

K - k LIM 

the correction maneuver will be initiated. A linear analysis of 
pending antiparallelism and parallelism was made to determine approxi- 
mate relations between K LBI , the CM3 momenta, and the commanded 
momentum. Appropriate relations are developed to establish the check 
factor limit for a given normalized total momentum requirement. 

An approximate method may be used to determine the impact of 
K on the orientation of the CMG unit momentum vectors and on the 
extent of CMG system saturation. As the saturation or the antiparallel 
condition is approached, the CMG momenta are driven to approach the 
commanded momentum in a predictable fashion since the vector sum of the 

commanded momentum and the CMG system momentum is maintained at zero by 
the ATM control law. 

The permissible resultant CM3 momentum vector orientations can 
be visualized by considering the projections of unit vectors along these 
momenta onto a plane normal to the commanded momentum vector, as is 


done m Figure 7 . The commanded momentum vector, which is coming out 
of the paper, is indicated by a small solid circle in the figure, and 
the intersections of the CM3 momentum unit vectors and the normal 
plane are denoted by larger circles. Dots and crosses are used to 
define CM3 momenta coming out of or going into the paper, and the 
projection of these CM3 momenta is described by directed, dashed lines 
Angles between the various momenta are designated by e,. 

Relations between the saturation ratio and the variable check 
factor K are developed for the cases of impending saturation and 
antiparallelism. To simplify the form and application of these 
equations, it is assumed that the angles e k between the comanded 
momentum and the CM3 momenta and the angles between the CMC momenta 
are sufficiently small to yield 

sin e k 1 «k 

€ k 

COS ev ~ 1 

x _ 2 

and terms of higher order than e| are neglected. 

Saturation 

The general CM3 orientation for the case of impending saturation, 
shown in Figure 7 , is characterized by an arbitrary triangle formed by 
the tips of the CM3 momentum vectors and pierced at its geometric 
center by the tip of the commanded momentum vector. The angles 
( e 4' e 5' € 6) between the commanded momentum and the CM3 momenta are 
related to the angles (e 1 ,€ 2 ,e 5 ) between the CM3 momenta by 



Figure 7* ~ Impending saturation. 



[ 2e l + 2 4 - e § 

1 ^ 

[ 2e ! + 2£ J - «fj 

> 

P e 3 + 2e l " e 2 


L _ 

J 


From Figure 7 and equations (6) and Cl?), the saturation 


(1 


io become: 



while the variable check factor is 


£ - H 


+ 6o + el 


Equations (15) and (l 6 ) give 


HT-i(8 + K) (1 

as the approximate relation between the variable check factor and the 
total CM2 system saturation ratio. Figure 8 presents the linear 
approximation for impending saturation (eq. (17)) plotted on the K 
versus HT parameter plot for a typical simulation case where anti- 
parallelism was avoided at K - 1.0 and HT = 1/3 and the control 
system proceeded to satiiration at K = l.o and HT = 1.0. 


Figure 8.- Linear approximation for saturation. 


case of impending anti- 


Antiparallelism 

The general CMG orientation for the 
parallelism, sketched in Figure 9, is characterized by a parallelogram 
formed by the end of the CMG momentum vectors and the conmanded 
momentum vector. The angle (e^) between the commanded momentum and 
the opposing CMG momentum is related to the angles between 

the CMG momenta by the condition 

s ? - i s 4 + 2 «§ - e£) (18; 

from Figure 9 . 



Figure 9.- Impending antiparallelism. 
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Tlie saturation ratio becomes 


HT - | - £(e§ + s| - eg) 

= $ - &*? - «?) 


while the variable check factor is given by 

K - 1 - |( £ i + e | + e §) 

- 1 ’ }(3«f + eg) 

Two cases are of interest and correspond to approaches from 
HT < 1/3 (antiparallelism) and from HT > 1/3 (parallelism). For 
the first of these, equation (18) indicates that for e h > e ± and 
elimination of e 1 from equations (19) and (20) yields 


HT ~ i (2 + K) + Jt 


if terms on the order of are neglected in equation (21), then 


HT ~ ^(2 + K) 


for HT < 1/3. For this case. Figure 10 presents this linear approxi- 
mation for impending antiparallelism (eq. (22)) plotted on the K 
versus HT parameter plot for a typical antiparallel approach case 
where spacecraft control was lost at HT = 1/3. The linear equation 
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CHAPTER V 


CORRECTION FOR ADVERSE CM3 ORIENTATIONS 
Correction Maneuver 

The basic approach to be studied in reducing K when K L3M 
is reached is to command an inner gimbal motor actuator to rotate a 
single CM3 spin vector to oppose a pending adverse condition and hence 
maintain control. The ATM control law automatically compensates for 
the change in the CM3 system momentum due to the inner gimbal rotation, 
and will thus bring the remaining two CM3 momenta out of the adverse 
orientation while continuing to provide control for the ATM cluster. 
Since the CMG system momentum is fed directly to the gimbal rate 
commands, the CM3 gimbal torquers automatically compensate for any 
change in this momentum which is not due to spacecraft momentum 
requirements (Appendix A). The net result is a rotation of the CM3 
momentum vectors away from the adverse orientation without the 
introduction of spacecraft attitude or rate errors. 

K-Factor Technique 

To correct for an adverse CM3 momentum distribution, the gimbal 
rates can be commanded to move one or more CM3 momentum vectors away 

from the adverse orientation. A deadband on K can serve to implement 
such a scheme. 

The simplest momentum reorientation command would use a single 
inner gimbal rate (say that for the largest inner gimbal angle, which 
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can be readily determined by comparison of the available cosine 
functions for the three inner gimbal angles) and would drive that 

inner gimbal angle toward its zero value, until an acceptable value 
of K is reached. 

Figure 11 illustrates a possible implementation of such a 
correction for adverse CM3 momentum distributions. The maximum inner 
gimbal angle and the check parameter K would be continuously or 
periodically determined. When K exceeds a preselected value, K L3M , 
the maximum inner gimbal angle is commanded toward its zero value at 

= ~K sin (PMAx) (27) 

until K is reduced by and a more favorable momentum distribution 
has been established. The correction maneuver is terminated when 

k < 

This CM3 reorientation also serves to reduce the largest CM3 
inner gimbal angle, and thus aids in improving the corresponding CMG 
effective gain. A known torque Hf) c sin 3max is applied to the space- 
craft during the correction. However, this torque is directly 
compensated for by the CM3 system inner loop and can be selected to 
be sufficiently small to cause negligible spacecraft errors. 
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Figure U. - Constant deadband K- factor correction. 
Isogonal Distribution 

The most effective control (Ref. 2) with the three available 
CMC momentum vectors is provided when these control momentum vectors 
have equal angles between themselves and are grouped in an umbrella 
configuration about the command momentum as illustrated in Figure 12, 
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Figure 12,- Isogonal distribution. 

The dot products of the unit momentum vectors are 

(W ’ ( h 2}) = ((^i) * (h 3 )j = ((hg) • (h 5 )j = cos 7 
so that the normalized angular momentum 

1112 • I + f[(( h i> • (»a) * (»i) • H * fc) • (h 5 »] 

becomes 

HT 2 = i + | cos 7 (50) 

while the check factor K from equation (ll) 

K = (( h l) ‘ ( h 2)) (( h l) * ( h 3))((h2) * (h 3 )J (31) 

become; 

K = cos-^7 (52) 

Eliminating cos 7 between equations (30) and (32) yields the relation 
between the isogonal check factor, KI, and the normalized momentum 
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Figure 13.- Variation of the isogonal check factor KI 
with the normalized momentum. 

Extension of K-Factor Correction to Isogonal Case 
A simple extension of the K correction technique to yield an 
isogonal correction appears possible by introducing a variable 



deadband for initiation of the largest inner gimbal angle rat, command. 
The isogonal correction satisfies the equation 

KI = i[3HT2 - if ( 3 M 

for all values of HI. Since HT? is available from the measured CM3 
system moment™ direction cosine vectors, the isogonal value for KI 
can be directly calculated from this expression. 

The resultant value of KI could serve as a lower bound for 
the K correction, which would now be implemented as shown in 
Figure Ik. A variable deadband defined by a lower limit KI and an 
upper limit where 


K ltm 3 KI + £K 


and ^ is a preselected constant, is introduced. When K exceeds 
K LIM> the maximum inner gimbal is commanded at 


Pi'o \ = -p_ sin s. 
V Pmay ) c 


This correction drives toward zero and continues until K 

reaches the isogonal value KI. 



Variable deadband K-factor 






CHAPTER VI 

APPLICATION SIMULATIONS TO AVOID 
ADVERSE ORIENTATIONS 

Both the constant and variable deadband versions of the K 
correction technique have been simulated for the ATM mission. A 
flexible ATM spacecraft, with seven bending modes incorporated 
dynamically and the remaining bending modes incorporated quasi- 
statically, is considered. Bending data and control law parameters 
and filters are taken as the ATM design values (Appendix C). The CMG 
servo loops are assumed to be ideal, and yield a CM3 transfer 
function of unity. Example cases are presented here to illustrate 
the spacecraft and CM3 system response for a constant 20 ft-lb 
disturbance moment representative of a worst case gravity gradient 
bias disturbance. Zero-momentum and orthogonal initial orientations 
for the CM3 system are considered. 

Constant Deadband K-Factor Correction 
Figure 15 shows the ATM and CM3 system response for a character- 
istic antiparallel approach starting from a zero-momentum orientation. 
The applied disturbance torque opposes CM3 1. Rate and attitude errors 
at the center of the ATM experiment package are given in Figure 15(a). 
The errors associated with the K correction are small in comparison 
with the steady- state errors produced by the applied disturbance 





Constant deadband antiparallel approach from zero-momentum 
orientation (20 ft -lb moment opposing CMC 1). 
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(c) Performance factors. 
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defined as the ratio of the control torque magnitude to the disturbance 


torque magnitude for 


GAIN = 

It) 


Gimbal angle time histories are presented in Figure 15(b). Note the 
inner gimbal maneuver of and the resulting rotation or ai of 
CM3 1 into a favorable position for control in Figure 15(b), and the 
small GAIN variations during the K correction in Figure 15(c). 
Additional performance indices for this case are also presented in 
Figure 15(c). The upper portion of the figure displays the variation 
of HT and K. It can be seen that HT remains approximately linear 
with time as expected for the constant applied moment. The correction 
parameter K increases until the deadband upper limit of K LIM =0.7 
is reached and then is driven away from its antiparallel value into a 
more favorable path to saturation. The lower portion of Figure 15(c) 
shows inner gimbal and control gain performance parameters. The inner 
gimbal performance parameter, P p , is defined by the relation 

p 3 = 1 - i(cos + cos | 3 2 + cos 3 3 ) av (38) 

and serves as a measure of the average inner gimbal angle excursion 
(Appendix A, eq. (A28). Its average value at saturation here is about 
0.160 and approximately corresponds to an average inner gimbal angle 
of 33 . Except for small transient variations initially and during 
the K correction, the GAIN factor remains at unity up to saturation. 
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Constant deadband saturation approach from orthogonal 
orientation (20 ft-lb moment alone: CMG l). 
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Figure 1 6.- Concluded. 



Ideal gimbal power consumption for the antiparallel approach is about 
5^1 watt-sec. 

Figure 16 illustrates another example for a constant disturbance 
torque along CM3 1 starting from an orthogonal orientation. The K 
correction now is not initiated until the CMS system is near saturation 
and then produces small control torque and spacecraft oscillations 
until saturation. In the actual AIM mission, desaturation would 
generally be carried out before the CM3 system is fully saturated and 
all control is lost. Hence, these oscillations near saturation may 
not occur. Ideal gimbal power for this example is about 218 watt- sec. 

Variable Deadband K-Factor Correction 
The extension of the K correction technique to incorporate a 
variable deadband with a lower bound defined by the isogonal value of 
KI, equation (?H could continuously minimize the CM3 system inner 
gimbal angles and could improve the CM3 system effective gain. 

Initial antiparallel approach cases for the constant 20 ft-lb 
disturbance moment were made to determine the effect of parameter 
variations in the correction technique on the spacecraft response and 
cm system perfoman-o factors. The correction parameters are jj c , 
the command gimbal rt e for the largest CMC gimbal angle and M 
which when added to the fixed isogonal value, KI, determines K LD1 
for initiation of the correction maneuver. For the antiparallel case, 
the spacecraft response was insensitive to changes in either B c or 
M as shown in Figure 1 7 . The gimbal power requirements, however. 
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increased with 3 C for all values of AK. Figure 18 illustrates the 
increasing power requirements and shows that 365 watt-sec of power was 

required for the characteristic antiparallel case with 0 C = 3 deg/sec 
and AK = 0.1. 
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The inner gimbal performance factor was also examined for 
correction parameter trade offs and is presented in Figure 19. 
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19*- Inner gimbal angle performance 
and AK with the variable deadband 


factor for variation of 
correction technique. 
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The data show for increasing command gimbal rate, p c , the average 
inner gimbal angle decreased for the antiparallel case. However, no 
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significant increase in the factor was noted in decreasing the 
command gimbal rate from 5 deg/ sec to 3 deg/sec. This fact, together 
with the trend of increasing gimbal power requirements with increase 
in 3 C , resulted in selection of an inner gimbal correction rate of 
Pc = 5 deg/sec as a standard to evaluate the overall performance of 
the K-factor correction technique. For the delta check factor deadband, 
a nominal value of £K = 0. 1 was selected and used with p c = 3 deg/sec 
for the following examples of the variable deadband scheme. 

Figure 20 again presents an antiparallel approach from the zero- 
momentum CM3 orientation; and Figure 23 illustrates the response for a 
constant moment applied along the minimum inertia axis starting from an 
orthogonal CM3 orientation. From Figures 15(c) and 20(c), it can be 
seen that the use of the variable deadband K-factor correction for the 
antiparallel approach improves the inner gimbal performance of the CM3 
system and significantly reduces gimbal power requirements. The ideal 
gimbal power consumption is now 365 watt- sec for the antiparallel 
approach. In addition, the variable K correction now no longer 
introduces excessive control torque oscillations near saturation. 

Figures 21 and 22 present the performance of this variable 
deadband scheme for the disturbance moment opposing CM3 2 and CM3 3, 
respectively. Both cases required the same ideal gimbal power require- 
ment of 365 watt-sec with the same average inner gimbal angle excursion 
of 28 as was the case for the disturbance moment opposing CM3 1. To 
permit the comparison of the K correction and the isogonal correction. 
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the antiparallel approach example was simulated with the isogonal 
correction. The results are shown in Figures 24, 25, and 26, and 
should be compared with Figures 20, 21, and 22 for the disturbance 
moment opposing CMS 1, 2, and 3, respectively. The spacecraft response 
is similar for the two cases, but the ideal gimbal power required by 
the CM! system increased by a factor of 2 for the isogonal corrections 
to approximately 795 watt-sec, while bhe inner gimbal angle performance 
index is reduced from 28° to 21°. Note that the control system main- 
tains an identical K versus HT distribution of momentum in addition 
to an identical average inner gimbal motion for all three isogonal 


cases. 






Figure 20.- Continued 
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Figure 20.- Concluded. 
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Variable deadband antiparallel approach from zero -momentum 
orientation (20 ft-lb moment opposing CMG 2). 
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Figure 22.- Continued. 
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Figure 22 . - Concluded. 











4 - - Isogonal antiparallel approach from zero -momentum 
orientation (20 ft-lb moment opposing CMO 1). 






Figure 24 • - Concluded. 




5.- Isogonai uriliparallel approach from zero-momentum 
orientation (20 ft -lb moment opposing CMS 2). 






6.- Isogonal antiparallel approach from zero -momentum 
orientation (20 ft -lb moment opposing CMG 3). 






Concluded. . 





CHAPTER VII 

APPLICATION SIMULATIONS TO 
EXIT FROM ADVERSE ORIENTATIONS 

of the CMC adverse orientations is the primary require- 
ment for optimal control with the CMC control system. In addition, 
however, an effective technique to avoid those undesirable CMS momentum 
orientations should also bring the system out of an adverse orientation 
if such an orientation is arbitrarily acquired. To examine the 
effectiveness of both the isogonal and simplified logic techniques in 
countering loss of control when the CMC system is placed in an adverse 
orientation, the cases of antiparallelism and anti saturation were 
considered. Initial gimbai angles and the command moment of 20 ft-lb 
were selected for these adverse orientations, aid then the performance 
of the two different techniques were evaluated in bringing the control 
system out of these orientations. 

Exit From AntiparaTLel Orientation 
With the control system positioned in an antiparallel orienta- 
tion, it was determined that the K-factor technique could successfully 
bring the system out of this adverse orientation and permit continued 
system operation. As Figure 27 (a) indicates, the spacecraft initial 
error is about 67O arc sec for the standard case of £K = 0.1 and 
0 C = 3*0 deg/ sec. Ideal power consumption up to saturation was 
!+86 watt- sec for this case and the average inner gimbai single was 3V 3 . 

For the K-factor technique, variation of M, or 0 C did not appreciably 
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Figure 28 . - Concluded. 
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affect the spacecraft response or gimbal power requirements for exit 
from the antiparallel orientation. 

The isogonal exit from the antiparallel orientation, presented 
in Figure 28, resulted in an initial spacecraft error of approximately 
4000 arc sec with large GAIN factor variations. In addition, the total 
gimbal power requirement was approximately 7 63 watt- sec, and 

the average inner gimbal angle of 31° was smaller than for the K-factor 
cases. 

Exit From Antisaturation Orientation 
To place the control system in the antisaturation condition, 
the initial conditions for the CM3 control system were chosen from the 
final gimbal angle positions given in Figure 16, and the sign of the 
command momentum was changed to M* = -20 ft-lb. 

Figure 29 presents data for recovery from the antisaturation 
orientation for the case of 3 C = 3 deg/sec and M = 0.1. The 
initial transient response results in an x-axis error of approximately 
480 arc sec with a GAIN factor of 4.9 before damping to a steady- 
state level of 140 arc sec. Under the continuous application of the 
bias moment, the control system again approaches saturation in about 
400 sec with a total power consumption of 510 watt-sec. For this 
K-factor case, the resulting average inner gimbal angle was 31°. 

Figure 3 0 presents data for the isogonal correction and indicates 
the spacecraft response and system performance is essentially the same 
as for the previous K-factor case. The isogonal correction provides 
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for a smoother transition to eventual saturation resulting in a total 
power requirement of k22 watt-sec with an average inner gim’cal angle 
of 28 °. 





Figure 29.- Continued. 
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29.- Concluded 










Continued. 





CHAPTER VIII 


CONCLUSIONS 

At ATM design gains the K-factor correction technique described 
in this section effectively avoids adverse CMS orientations up to 
saturation and produces control for exit from the adverse orientations 
in the advent they are arbitrarily acquired. The correction thus 
improves the CMS system effective gain and, by reducing the CMC inner 
gimbal angles, requires less power than the H-vector control law alone. 
The K-factor correction can be readily implemented by analog or digital 
computer, and is an order of magnitude less complex than the current 
baseline isogonal correction. For the 20 ft-lb disturbance moment 
example computer simulations, the K-factor correction technique 
required considerably less CMG gimbal power than the isogonal correc- 
tion, with the exception of exit from the antisaturation orientation, 
but had slightly larger inner gimbal angle excursions. 
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APPENDIX A 


CONTROL SYSTEM EQUATIONS 


Control Law Development 


The control system for the AIM spacecraft uses a modified transpose 
steering matrix and CM3 system momentum feedback (Ref. A). The basIc 


governing equations for this system are reformulated here to Indicate 


the assumptions made In simulating Am control with the antiparallel 


correction method. 


Hardware components, such as the CMS units and the spacecraft 


rate and attitude sensors, are assumed to have a transfer function of 


unity and thus become Ideal components during the computer evaluation 


of the antiparallel correction. Torque contributions due to glmbal 


inertia terms are also neglected. The CM! system control torque {g} 
can thus be expressed in terns of the CMG momentum vector (¥} and the 
ATM spacecraft angular rate vector {ji) by 


(0> . -<H> - H(h) 


where 


(g) = /g 2 \, (h} = /H 2 l (ft) = Jn. 


R 
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and the cross-product term ^ 0 ^ is given by 

0 ^2 

Ql] = ^3 0 -f2^ 


■-02 fip 0 


In index notation, the ith component of the control torque form 
equation (Al) is 

Gf = H((aj_ •* ^)sin cos 3j_ + (a^ - Hj)cos cos 0^ 

+ (3i cos + Gj)sin 3 t - (3 k sin + G k )sin 3 k 
+ (3^ + flj sin aj + ^ cos a.j)cos 3j) (a* 

in its expanded form. Since the product of vehicle rates and CMG 
momenta is small, equation (Al) can be approximated by 

<°> “ ■{»} (a; 


for the purposes of control law analysis. 

The CMG momentum vector in equation (A5) is 

cos 3q cos -cos 3j sin 

{h} = H< cos 32 cos ®2 -cos 3q sin 

cos 3^ cos -cos 35 sin Or, 

so that equation (A5) can be rewritten as 


-sin 3p 
-sin 33 
-sin 3^ 


( a ) = -h(a] ( 5 ) 


where 


with 



MM 



sin 3j_ cos 
-sin (3 -l sin 0]_ 
cos Pj 


cos ^2 

sin p2 sin &2 
-sin pg sin 02 


-sin Pj sin aj 
cos Pj 
sin Pj cos 



sin 03. 0 

< 

] ■ Cm] Cm] - 

cos a]_ sin 02 

C 


a 3 


COS P]_ 
0 


0 0 
cos P2 0 
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The gimbal rate vector 6 for the ATM system is commanded to be 



or in index notation 

a-i = K SIj(( H Ci “ %)sin + (Hcj - Hj)cos 04, 

0i = %L^®Ci “ %)cos sin Pj, - (H^j - Hj)sin sin p^^ (AU-) 

- (Hck " H k) cos 0k) 

with the result 

(g) = -Ks L H [&d MF + M 00 M T J ({He) - (H>) 

* -% L h|b]((Hc> - (H>) (A12) 

The two 3*3 matrices making up [V] in equation (A12) are reasonably 
diagonal and permit the generation of continued control commands in the 
event of failure of a CMS unit. The gain factor Ksl is a constant. 

The commanded momentum He in equations (A10) and (A12) is 
derived from desired control torques of the form 



where (fi) and (q) are developed in Appendix B. The gains Ki are 
constant. A commanded moment is obtained from {mc} by fourth-order 
filters, which eliminate high frequency input to the CMG system. The 
filtered commanded moment is thus 



(A14) 


where s denotes the Laplace transform variable. 


The commanded momentum 


becomes 
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The magnitude of the commanded moment in equation (A15) is ^ 

limited to ±200 ft-lb because of the integrator limitations in the ATM 

control computer. , 

I Equations (A5) and (A12) yield 

(h) = K sl h[b] ({He} - (H)j (Al6) 

|g or 

l 

i (h) + k sl h[b] (h) = k sl h[b] (He) (ai 7) 

■: For lar e e values of K SI H, the effective response of the system 

described by equation (A17) is given by 

j. ( H ) {** 0 } (ai8) ' 

and 

.4 ( G ) *{%} (A19) < 

k 4 

Since the CMG system momentum is fed back to the CMG gimbal i 

torquers, any change in this momentum without a corresponding change 

in the commanded momentum must be compensated for by a redistribution 

of the CMG unit momentum vectors to maintain the condition of 

equation (Al8). This inner-loop response of the ATM control law is 

extremely fast (with a time constant of l/K SL H) and should be able to j 

compensate for such a CMG system momentum change before the effects 

# 4 

>; of the change are transmitted to the spacecraft. 

i 

i 
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Inner Gimbal Roll-Off Logic 

Also incorporated into the control law logic was a roll-off 
acheme for the inner gimbals when they reached their ±8u° stops. When 
an inner gimbal saturates as such, the corresponding outer gimbal is 
commanded to move in a direction as to back the inner gimbal off its 
stop to a specified liait (selected as p iL = ±75°) and at a specified 
outer gimbal rate (selected as ai c =5-5 deg/sec). 

Performance Factors 

Several factors are used to permit assessment of the CMG system 
performance. The first performance factor is the theoretical gimbal 
power consumption P . The electrical power Pqi for tne ith gimbal 
torquer is approximately expressed by the relation 

P Gi = ^B^i + i2r (A20) 

where tie first term represents mechanical power and the second term 
electrical power losses. Here 5± is the gimbal race, K B reflects 
the back EMF constant, Kq is the gear ratio, and R is the winding 
resistance. The current I is given by 

T i 

1 = Kt*g ( A21 ) 

and is a function of tne gear ratio Kq, the stall torque constant Bp, 
and the effective gimbal torque command Tp. Substitution of 
equation (A2l) into equation (A20) yields 
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P “ = + (A22) 

The effective control commands requiring power for the outer and inner 
gimbal become 


= (HPi cos + Tp sgn ai) 
T P i = cos 3 i + Tp sgn 3^) 


with 


(A23) 


T ai - 0.0 if sgn T ai ± sgn a^ 
T 3 i = °*° lf sgn T p _, £ sgn 0. 


(A24) 


where Tp is the gimbal friction torque magnitude and gimbal 
acceleration torques are neglected in equation (A23). When T ^ is 
negative, existing precessional torques are already producing tne 

desired control moments and the effective torque command for that 
gimbal goes to zero. 

Equations (A22) and (A23) are used to define the total CMG system 
gimbal power as 
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and the power consumption becomes 


P 



P G dt, 


watt-sec 


(A26) 


to complete the definition of the first performance factor. 

The second performance factor GAIN is defined as the ratio of 
the control torque magnitude to the disturbance torque magnitude, so 
that 


GAIN = — 
M D 


(A27) 


This factor ser/es as a measure of the effective gain of the control 
system. 

The third performance factor provides an indicator of the extent 
of saturation of the inner gimbals, and is defined as 

p e “ 1 ' 3 Z cos p i 
i =1 


(A28) 


APPENDIX B 


SPACECRAFT EQUATIONS 


Rigid-Body Contributions 

The rigid-body rates for the ATM spacecraft are determined 
integration of* the Euler equations 


(%} - QJ 1 ^) ♦ <%} - [%][i]{%>) 


where (g} is the control torque defined by equation (Al), {m d } is 
the applied disturbance moment, and the rigid spacecraft inertia and 


angular rate are defined as 


0- o 


0 0 


0 0 


° %1 
0 ' {%} H < n R2 


Products of inertia are neglected. 

The rigid-body Euler rates become 


Or) = \®R2/ * 
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Vhere 


1 



0 


tan 0^2 sin 
cos 

sec 0^2 sin e R1 


tan 0^2 cos 0j^ 
-sin 0 R1 
sec 6^2 cos 


(Bk) 


gives the transformation matrix between spacecraft and Euler rates. 

Integration of equations (Bl) and (b 4 ) defines the rigid-body 
motion of the ATM spacecraft. 


Flexible-Body Contributions 

The flexible -oody rates and attitudes are determined from 
dynamic and quasi-static modal representations of the ATM cluster. 
Seven modes are represented dynamically; all remaining modes are 
incorporated quasi-statically. The dynamic modes contribute both 
angular rate and Euler angle errors, while the quasi-static modes 
contribute only Euler angle errors. 

The basic uncoupled modal differential equations for the dynamic 
modes are 

<%>} + - &>]{ E ) (B5) 

or in index notation (for i =1,2, . . ., 7) 

.18 

<*Di + = Z U Dij E j (Bo) 

J=1 


where q Di is the modal displacement, % is the modal damping, and 
“bi is the modal frequency. The vector (e) contains the l8 moment 
and force inputs (defined in Appendix C) corresponding to disturbances 
applied in the ATM/lM (Ep to Eg), the CSM (Ej to E 12 ) , and the SIVB 
Workshop (Epj to Epg) . The matrix (^Ud^ is the 7 x l8 modal 
distribution matrix for the seven dynamic modes considered. The 
flexible raoes and attitude errors for the dynamic modes become 

{%} = [0d]( 4 d > (® 7 ) 


( 7 fd) " M^ D ) 


or in index notation 


o o 

r FDi = Z ^Di^i. = Z U Dji^i 


where q$ and q D are determined by integrating equation (B5), and 

M 3x7 modal participation matrix found by transposing the 

first three columns of [U;q] . 

For all remaining modes, which are represented quasi-statically^ 
the modal displacements become 
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by neglecting terms containing derivatives of the modal displacement ii 
the modal differential equations , the resulting attitude errors are 


( r Fs) " 


5 DOC*} 


The total flexible-body attitude errors correspond to 


{ 7 f) { 7 fd} + { 7 fs} ( B1 

from equations (b 8) and (Bll) and the flexible Euler angle errors are 


given by 


( e s) - Q0<7 F } (B13) 

where QQ is determined in equation (B^). 

Evaluation of equations (B7) and (Blj) defines the flexible -body 


motion of the ATM spacecraft . 


Total Spacecraft Errors 


The total spacecraft errors are found by adding the rigid-body 
and flexible-body contributions to the spacecraft motion. This gives 


(n) - + {n F ) 


using equations (Bl) and (B7) for the rate 


errors and 


(e) = (e R ) + (p f ) 


(B15) 


using equations (b 4) ana (B13) for the Euler an«l 


angle errors. 




APPENDIX C 


SIMULATION CONSTANTS 


The program constants used in the ATM computer simulations are 
included in the following sections. For the ATM cluster the moments of 

inertia were chosen as: 

lx = 595151 s.’.ug-ft 2 

Iy = 5866985 

I z = 5755771 

with zero products of inertia. The integration computing interval was 
0.05125 sec for the equations of motion and the control law command 
signals. The ATM design gains were taken to he 


Kp = 29500. ft -lb 

Kj =» 11+7512. 

K5 = 568781. 


K 2 = 16226k. ft -lb -sec 

K4 = I505I+85. 

Kg = 18I+3906. 


%l = 0.00525 per ft -lb -sec' 


In the control system equations the constant terms for the 
fourth order filters were 


P 1 * p 2 “ P3 “ 0.5 
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for the damping factor and 


= 1.25 rad/sec 
co 2 = 1.25 
o > 5 = 1.90 

for the filter frequencies. 

flexible spacecraft modal differential equations include 
the modal frequencies 

“U! = O.9362 

U>Q2 = 1.0116 

Uijjj = 2.6012 
= 3.3929 
u>)5 = 4.0087 
o>OS - 13.6534 

Ufff = 15.8336 

and the modal damping 

£ = 0.01 


with the forcing function 
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> 


C l8 J 

modified by the modal distribution matrix QjJ given in Table I. 

The quasi-static contribution to the flexible spacecraft response 
requires the residual compliance matrix [V] given in Table II. For 
the flexible spacecraft attitude solution the modal participation matrix 
[f] is given in Table III. 


In the gimbal power equations the following values were used 
for the ATM torque constants. 


Gear ratio = 56 

Back EMF Kg = 1.53 volt/rad/sec 

Stall torque % = 1.15 ft-lb/amp 

Winding resistance R = 6.9 ohms 
Friction torque Tp = 5.0 ft-lbs 
For the isogonal correction (Ref. 2) a distribution law gain of 

Kj) * 0.1 and a rotation law gain of Kr »-0.01 were used for the AIM 
mission simulations. 
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